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. We estimate the Cronin effect in pA collisions at the CERN LHC and at RHIC, using a 
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Glauber-Eikonal model of initial state multiparton interactions. For a correct determination 
of the initial parton flux, we upgrade the model cross section, taking carefully into account 



1. INTRODUCTION 



all kinematical constraints of each multi-parton interaction process. As compared with pre- 
vious results, derived with approximate kinematics, we obtain a softer spectrum of produced 
partons, while improving the agreement of the model with the recent measurements of 7r° 
\ production in d + Au collisions at y/s = 200 AGeV. 
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Hadron-nucleus interactions represent an intermediate regime between hadron-hadron and 

nucleus- nucleus collisions, where different ideas on the role played by the complexity of the 

target on large pt dynamics may be tested [1]. Of particular interest is the study of the Cronin 

^ effect [2, 3], namely the deformation of the transverse spectrum of particles, or jets, produced 

H ' 

in hadron-nucleus collisions in the projectile fragmentation region. Given the large scale of 
momentum transfer in the process, the problem is within reach of a perturbative approach, while 
the effects of the complex structure of the target may be controlled by changing energy and 
atomic mass number. 

The effect is originated by the high density of the nuclear target, which induces the projectile's 
partons interacting simultaneously with several partons of the target. There is not however a 
unique explanation of the underlying dynamical mechanism [4-8] and different approaches can 
only be discriminated by measuring the large pt spectra in hadron nucleus collisions at the 
energies of the future colliders, where various models give different predictions. In our opinion, 
the simplest approach is to adopt the Glauber prescription of factorization of the overall many- 
parton S matrix, as a product of elementary partonic S matrices [9]. Since production at the 
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parton level represents a higher order correction in the coupling constant, the natural zeroth 
order approximation is obtained when evaluating each elementary parton interaction at the 
lowest order in as, which leads considering binary processes only. 

Encouragingly a recent calculation, based on this simplest scheme [10], compares rather rea- 
sonably with the latest experimental measurements of the Cronin effect in large pt 7r° production 
in d + Au collisions, at a center of mass energy of 200 AGeV [11], [12], [13], [14], which represents 
a convenient channel to study the effect (as, differently from charged particle production, there 
is no contamination from protons, whose production mechanism poses additional problems). 

It's remarkable that in such an approach, when adopting the usual approximations of the 
Glauber picture of high energy hadron-nucleus collisions, the whole effect of the nuclear target 
medium may be taken into account exactly, in the evaluation of the inclusive transverse spec- 
trum [9]. After summing all possible multi-parton interactions of the projectile parton one in 
fact obtains 

where 



e T A (b)[F A (x,r)-F A (x,0)} _ e ~T A (b)F A {xfl) 

and 2~k(6) is the usual nuclear thickness, as a function of the hadron-nucleus impact parameter 
b, G(x) the parton number density of the projectile, as a function of the fractional momentum 
x, pt the transverse momentum of the final observed parton and 

FA(x > b) = I i^ dx ' GA{x,) w/ ptb (2) 

In Eq.(2) Ga(x') is the parton number density of the target nucleus, divided by the atomic 
mass number A, as a function of the fractional momentum x' , and da/d 2 pt the elementary parton- 
parton scattering cross section, which includes also the kinematical limits constraints. For sake 
of simplicity flavor dependence is not written explicitly. The quantity S^ ard (r, b, po) can be 
interpreted in terms of dipole-nucleus hard cross section. The dipole originates from the square 
of the scattering amplitude, expressed as a function of the Fourier variable f , corresponding to 
the transverse size of the dipole [15]. 

By expanding Eq.(l) as a power series in a, one obtains all different multiple scatterings of the 
projectile parton. As a consequence of the Glauber prescription of factorization of the overall 
many-parton S matrix, each term in the expansion is represented by an incoherent convolution 



S hard( r i b > PO) = 
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of binary partonic collisions [9, 16]. Analogously to the case of the Glauber approach to hadron- 
nucleus processes, the 5-matrix factorization prescription does not however imply any space- 
time ordering between the different elementary collisions. It only implies that a connected 
ra-body interaction process is well approximated by a product of two-body interactions, which 
corresponds to the dominance of the pole contribution, in a dispersive representation of the 
projectile-exchanged gluon amplitude [17]. 

All expressions are divergent for small values of the lower cutoff in pt, which sets the regime 
of hard interactions. Nevertheless, as unitarity is explicitly implemented, each term in the series 
expansion of Eq.(l) contains subtractions, which lower the degree of the infrared singularity, from 
an inverse power to a power of a logarithm of the cutoff [16]. It's also worth mentioning that, 
in this model, the average energy, carried by the partons scattered to the final state is bound to 
be smaller than the overall initial state energy [18], whereas such a quantity is divergent as an 
inverse power of the cutoff, if evaluated in impulse approximation. 

The approach to hard processes in hadron-nucleus collisions, based on the factorization of the 
multi-parton S matrix, is hence unitary and satisfies various non-trivial consistency requirements, 
including the AGK cutting rules [19]. As pt — > the importance of unitarity corrections grows, 
producing a suppression of the integrated parton yield and a random walk of partons to higher 
Pt , recovering in this way the local isotropy in transverse space of the black disk limit, which is, 
on the contrary, maximally broken by the lowest order impulse approximation term [18]. 

The model hence contains geometrical shadowing at low pt and Cronin enhancement of the 
transverse spectrum at moderate pt- Some non secondary features of the interaction are not 
accounted however, as one expects genuine dynamical shadowing at low x, due to non-linear 
gluon interactions [1]. In spite of some interesting efforts in this direction [20], incorporating 
non-linear gluon evolution in the present formalism is however a non trivial problem, which we 
will not try to address in the present note. 

In fact we aim to point out that also some of the simplifying assumptions, adopted in the 
derivation of Eq.(l), require a closer inspection. More specifically, to obtain Eq.(l), one needs 
to decouple longitudinal and transverse degrees of freedom, while conserving the longitudinal 
momentum components in the interaction. Kinematical constraints are hence implemented only 
approximately, in such a way that final state partons turn out to be more energetic than the 
initial ones, the effect being emphasized when the number of re-scatterings grows. The final pt 
spectrum gets hence shifted towards larger transverse momenta which, given the steepness of 
parton distributions at small x, might modify appreciably the evaluation of the Cronin effect. 

The purpose of the present note is to investigate this specific aspect of the model more closely. 
We hence compare the spectrum, obtained by implementing exactly all kinematical constraints, 
with the summed expression in Eq.(l), both at the LHC regime and in the kinematical domain 
of the recent measurements of the Cronin effect in d + Au — ► ir°X at RHIC. 
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2. TRANSVERSE SPECTRUM 

The expression in Eq.(l) represents the sum of all possible collisions of the observed parton 
with the different scattering centers of the target nucleus. When the kinematical constraints 
are implemented exactly, the multi-scattering series cannot be re-summed any more in a closed 
analytic expression and each different term needs to be evaluated separately. Most of the spec- 
trum is nevertheless well reproduced by the the first three terms of the expansion [16]. Hence 
we approximate the cross section by: 



da s-^ doj 

J2n*d„J2h { > 



d 2 ptdyd 2 b ^— ' d 2 ptdy d 2 b 

where pt is the transverse momentum of the observed parton, y its rapidity, b the impact param- 
eter of the collision and the sum runs over the different species of projectile partons, while 



da, da[ 1] da[ 2) da^ 



j 1 1 1 (4) 



d 2 pt dy d 2 b d 2 pt dy d 2 b d 2 p t dy d 2 b d 2 pt dy d 2 b 



where 



, (i) 1 /" 

x x \ fi (x, Qf ct ) f ± (x', Q fct ) + f ± (x 1 , Qf ct ) fj_ (x, Q fct ) 

Lp A p A 

is the single scattering term (the projectile parton interacts with a single parton of the target 
and vice- versa), 



d?p t dyd 2 b ... , r u n 

x 5 (2) (qi +q 2 - Pt) - 5 (2) (gi - Pt) - 5 {2) (g 2 - Pt) 

x xfi(x, Qfct)fA x 'i Qfct)fk(x", Qf ct )afl{y,x' \x"\qi,q 2 ) 

p A A J 

is the rescattering term (the projectile parton interacts with two different partons of the target), 
and 



5 



= a YTs~ TTs~ k TT5~i J ^ d ^ d ^ 

3 ,k,l J 
x dx dx" dx" af) H (y,x',x", x" ; qi , q 2 , q 3 ) 

x <5 (2) (qi + q 2 + qz~ Pt) - 5 (2) (qi + q~2 - Pt) - <5 (2) (q~2 + <?3 - Pt) 
-5 (2) (? 3 + Qi ~ Pt) + 5 {2) (qi ~ Pt) + 5 (2) (02 ~ Pt) + 5 (2) (g 3 - ft) 

X Qfct)f±(x', Qfct)f±(x", Qfct)f±(x"', Qfd). 

p A A A 

is the double rescattering term (the projectile parton interacts with three different partons of 
the target). 

The subtraction terms, in the rescattering and in the double rescattering terms, are a direct 
consequence of the implementation of unitarity in the process. The sums run over the different 
species of target partons. The quantities f±(x', Qf c t) represent the nuclear parton distributions, 

A 

while the multiparton interaction cross sections are given by: 

a^(y,x'; qi ) = k factor ^ 8(*W) 

dU 



J dto dU 



^ f„ rrf n» ^".n n n \ - (h ^3 d ^il( § 3, k, U 3 ) dd ik (s 2 , t 2 , U 2 ) 
a iikl\y> X > X > X ; 51)92, 03) - {^factor) ^ ^ 

J dfa dt 



3 az 2 

e($(3)) ($(2)) ( $ t 



daijjh, ii, Ul) Bf <|)( ; 5) , H( (T) (2i , f _ V(|)( J , - 



where the partonj-parton., differential cross sections d&ij/dt are evaluated at the lowest order in 
pQCD by making use of exact kinematics. The kinematical constraints are implemented through 
the quantities . The expressions of $W and of the invariants lj , ii, Ui are derived in great 
detail in the appendix. 

The infrared singularity is regularized by introducing in the propagators the mass parameter 
Po , which is a free parameter in the model. Higher order-effects in the elementary interaction 
are accounted by multiplying the lowest order expressions in a s by the factor fc/ ac t- 

When evaluating the 7r° spectrum, we let partons fragment independently, neglecting the 
transverse momentum component generated by the fragmentation process. The relations between 
hadronic and partonic variables are hence the following [21]: 
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m h cosh y h = zp t cosh y m h sinh y h = q t sinh y 

where z = ^ , with £7 the energy of the parton, y its rapidity, Eh the energy of the hadron h, m 
its mass, yh its rapidity, qt its transverse momentum and mh its transverse mass. The inclusive 
7T° spectrum is hence given by: 

da^ _ cfcr/ rafl 
d 2 q t dy h d 2 b ^ d 2 q t dy h d 2 b 

with 



daf ra9 , f dz _ , „o , do; 



/dz 



d 2 q t dyhd 2 b J z 2 ' d 2 p t dyd 2 b 
where 

Qt T / \ ( m h sinhw 

Pt = — J {m h , yh) V = arcsin ' 



and Di^h{z,Q 2 F ) the fragmentation functions, depending on the fragmentation scale Qp. The 
integration region for the hadron energy fraction z is set by the request pt > po, with po the 
infrared regulator: 



2m h coshy h 



< z < min 



1, — J{m h ,y h ) 
Po 



2.1. Numerical Results 

In our calculations we use the isospin averaged nuclear parton distributions, f± = Z f± + 

A p 

(A — Z) f± , where for the proton and the neutron distributions, /; and f±, we take the CTEQ5 

n p n 

parametrization at the leading order [22], while for the nuclear thickness function Tk(6), which 
is normalized to one, we use the Wood-Saxon parametrization. 

To study the effect at the LHC we consider the case of production of minijets in a forward 
calorimeter (r? 6 [2.4, 4]) at two different center of mass energies in the hadron-nucleon cm. 
system, = 5.5 ATeV and yfs = 8.8 ATeV, which might be compatible with the acceptances 
of ATLAS and CMS [23]. Given the small values of x in this kinematical range, we use effective 
parton distribution functions, with the gluon-gluon cross-section as a dynamical input. We set 
the factorization and renormalization scale Qf c t, Qrn equal to the regularized transverse mass 
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m t = y/l>o + Vt an( i we evaluate the transverse spectrum of jets, in proton-lead collisions, by the 
expression: 



dW ' A 1 f P) 



(Ppt Vmax-Vmin J yd[y mln ,y max ] ' d?p t dy d?V 

where the integrand is given by Eq.(3). As for the choice of the infrared regulator we set 
Po = 2GeV and we use kf ac t or = 2. Our results are plotted in Fig.l, where we compare 
the spectrum (5) obtained by the exact implementation of energy conservation in the multiple 
interactions (solid line), with the approximate kinematics results given by the first three terms 
of the expansion of Eq.(l). As an effect of the exact implementation of kinematics the spectrum 
of the outgoing particle is shifted towards lower transverse momenta, resulting in an appreciable 
energy loss effect: at pt ~ 15 GeV the energy-loss correction to the spectrum is of the order of 40- 
50%, at higher transverse momenta, pt ~ 30 GeV, it's still about 30-38%. The triple scattering 
approximation, used to evaluate the spectrum, breaks down at pt < 9 GeV at y^i = 8.8 TeV: 
by increasing the center of mass energy the density of target partons grows rapidly and the 
contribution of higher order rescatterings cannot be neglected any more at pt < 9 GeV (right 
panel of Fig.l). As the energy is lowered to yfs = 5.5 TeV , the region of numerical instability is 
shifted to the region pt < SGeV (left panel of Fig.l). 

The softening of the spectrum has important effects on the Cronin ratio, defined as: 

dW pA /d 2 Pt 



R(pt) 



dwl l }/d?p t 



where dW^ /d 2 p t is the cross section in impulse approximation. The Cronin ratio, corresponding 
to the spectra in fig.l, is shown in fig. 2, where the dashed and continuous lines correspond to 
approximate and exact kinematics respectively. The two panels of fig. 2 refer to the two values 
considered for the center of mass energy; at higher energy multiple interactions produce a larger 
value of the Cronin ratio. The overall effect of the implementation of energy conservation on 
R is to decrease the size of the Cronin effect and to shift the curve towards smaller transverse 
momenta. 

At lower energy, yfs = 200 AGeV, we follow [10] in the evaluation of the cross section of 
d + Au — ► tt°X, using Qf c t = Qrn = ^ , and the values po = 1.0 GeV and fc/ ac t = 1-04. These 
choices, besides minimizes the effects of higher order corrections in as, allow one reproducing 
the inclusive cross section of 7r° production in pp collisions at the same cm. energy, without any 
need of smearing the cross section with an intrinsic parton transverse momentum. At ^fs = 200 
AGeV we take explicitly into account the differences between partons, both in the distributions 
and in the evaluation of the partonic cross sections. The resulting Cronin ratio in d + Au 
collisions is hence a parameter free prediction of the our model. By using the leading order 
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K-K-P fragmentation functions [24] at the fragmentation scale Qf = 



we evaluate: 



RdAu-*-K°X — 




I 




d 2 qt dy d 2 b 



cPqt dy d 2 b 



at y = and b = b^Au = 5.7 fm, which is the estimated average impact parameter of the exper- 
iment [25]. In fig. 3 we compare our result (continuous line) with the experimental data of the 
PHENIX collaboration [11] and with the result obtained by using the approximate kinematics of 
Eq.(l), dashed line. Given the small rapidity values of the observed 7r°, in this case the correc- 
tions induced by exact kinematics are particularly important in the region of smaller transverse 
momenta. The dependence of the effect on the impact parameter is shown in fig. 4, where our cal- 
culation (continuous line) is compared with preliminary data presented by PHENIX at the DNP 
Fall Meeting, held in Tucson last October 2003 [26], and with the the standard Glauber- Eikonal 
calculation (dashed line). The fractions refer to the BBC centrality selection. The corresponding 
estimated impact parameter values are: 0-20% => b=3.5 fm 20-40% — > b=4.5 fm 40-60% — ► 
b=5.5 fm 60-88% — > b=6.5 fm. As an effect of implementing exactly kinematical constraints, 
one obtains a reduction of the Cronin curve in the region of smaller pt values, resulting with 
an improved agreement of the model with the experimental indication. Notice that besides the 
deformation of the spectrum leading to geometrical shadowing, by implementing energy conserva- 
tion one obtains also the quenching the spectrum, as part of the projectile's energy is transferred 
to the target. 

The expectation for tt° production in d + Au at rj = 3.2 is shown in fig. 5, where we compare 
our result (continuous line) with the standard Glauber- Eikonal model calculation (dashed line). 
Given the larger average number of rescatterings at larger rapidity values, the quenching of the 
spectrum, due to the energy lost by the projectile in the multiple collisions, is now sizably in- 
creased (notice that the vertical scale is now different with respect to the previous figures 3 and 
4). One in fact obtains values of RdAu-*it°x as small as .6 at pt = 1 GeV. The expectation is 
nevertheless that the Cronin curve should exceed one for p t > 2 GeV, which is not the case for 
charged particle production. In the figure we plot preliminary results for charged particle pro- 
duction at 7] = 3.2, presented by BRAHMS at the Quark Matter Conference[27]. In the forward 
region the quenching of the spectrum might hence be a sizably larger effect, not accountable by 
multiple scatterings. 



In very high energy proton- nucleus collisions the interaction probability is very large, also 
when considering processes with momentum transfer well above the lower limit for a perturbative 
approach. In this regime unitarity needs to be implemented already at the level of perturbative 
interactions and the Glauber-Eikonal model of multiple parton collisions allows implementing 



2.2. 



Conclusions 



9 



unitarity in the simplest and most natural way. A standard approximation is to conserve the 
longitudinal momentum component of the projectile in the process. The approximation becomes 
however questionable when dealing with parton interactions, since parton distributions are sin- 
gular at small x and a small change of the parton momentum leads to sizable modifications of 
the initial parton flux. 

In the present note we have studied this particular aspect of the problem, working out the 
leading terms in the expansion in multiple parton collisions and taking into account all kine- 
matical constraints exactly. Our result is that the spectrum of the produced large pt partons is 
appreciably softened after implementing exact kinematics. 

The Cronin effect has been recently measured in the large pt spectrum of 7r° production in 
proton-lead collisions, at y = and \fs = 200 AGeV. We have worked out the 7r° spectrum in the 
same kinematical conditions and found a rather good agreement of our model with experimental 
data. In this case the corrections to the Glauber Eikonal model, induced by implementing 
kinematics exactly, are increasingly important in the region of smaller p t , where, given the 
rapidity value of the observed 7r°, longitudinal and transverse momentum components are of 
comparable magnitude. The quenching of the spectrum, due to the energy lost by the projectile 
through multiple collisions, is a stronger effect at larger rapidity. We have estimated the effect 
at i] = 3.2, where the preliminary experimental indications seem however to indicate a stronger 
quenching effect than expected in the model, which, at relatively large pt, still gives a value of 
the Cronin ratio larger that one. 
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4. APPENDIX 

The four-body parton interaction (one projectile parton against three target partons) is rep- 
resented in the model by the convolution of three on-shell elastic parton-parton scatterings. We 
call h the four-momentum of the initial state projectile parton, and p' , p" ,p"' , the four-momenta 
of the three nuclear target partons. All partons are massless and on shell and have zero trans- 
verse momentum components. By introducing the momentum fractions x, x' , x" , the light-cone 
components of the four momenta in the initial state are: 
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h= (Vsx,0;0) p' = (0,y/sx';0) p" = (0,x/sx";0) p'" = (0, y/sx'"; 0) 

The transverse components of the exchanged four-momenta are ft, (72 and 53, while /, J, if 
are the projectile's four-momenta after the first, second and third interaction and P, Q and R 
the four-momenta of the recoiling target partons. For the transverse components one has: 



I± =qi J± = qi + Q2 K± = ft + q 2 + 93 Pi. = ~Ql Qa 



-Q2 R± 



-Q3 



While the energy conservation and mass-shell conditions, for each two-body scattering process, 
are expressed by the following equations: 



I+I- 
P +p- 

s/sx 

y/sX 1 



ft 

I++P+ 
I-+P- 



(6) 



J+J- 
Q + Q- 
1+ 

I' + yfsx" 



(Qi + Q.2Y 



?2 



J + + Q A 
J- +Q- 



(7) 



K+ K- 
R + R~ 
J+ 

J- + yfsx'" 



(qi + Q2 + qzY 

K+ + R+ 
K~ + R- 



(8) 



The three systems of equations allow to determine 

x , I+,r , J + , J-,K~,P + , P- Q+, Q-,R + , R~ 

as a function of ft, q 2 , Q3, K + , x' , x" and x'" . 

From the system (8) and the mass shell condition J + J~ = (ft + q 2 ) 2 , one obtains: 



(y-sx"> K+ - (ft + q 2 + q 3 f) (J + ) 2 - (ft + q 2 f (K+f 

- (yfsx'" K + - 2 (91 + q 2 ) • (91 + <12 + &)) J + K + = 
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which has real solutions when A^ 3 ) > 0, where 

A( 3 ) = (^x"'K + ) 2 -A(q 1+ q 2 ).q 3 ^x'"K + 
+4 (((91 + g 2 ) ■ Q3? ~ (?i + 9 2 ) 2 qfj 

The condition A^ 3 -* > can be satisfied only if <E>( 3 ) > 0, where: 



$ (3) = i_ 2 ((gi +92) • 93+ I gi +g2 | | 93 

y/8x"'K+ 

which gives the kinematical limits for x'" and for K + 



ill = 2 ((91 + 92) -93+191 + 921193 |) = jC|n . 



y/8X"'K+ K+ 

One hence obtains J and R: 



J+ _ K+ x'" - 2 ((ft + 92) • (9x + 92 + 93)) + 



2 

in 



r = 
r + = 

BT = 



2 ^^ + ^"-fe+?2 + g 3 ) 

VaW + 2 ((9! + 92) • (91 + g 2 + g 3 )) - V~sK+ x 
2K+ 

K + y / A(3T + 2 (q 3 • (qi + q 2 + q 3 )) - yfi K+ : 

2 K+ x'" - (q! + q 2 + q 3 ) 2 

VaW + x'" - 2 (9-3 • (91 + 9"2 + 93)) 



2K+ 

while the Mandelstam invariants S3 = (J + p"') 2 ,i 3 = {K — J) 2 , u 3 = (R — J) 2 are expressed 
by: 



S3 = Jlx'"j + 

, _ ((gi + 92 + 93) J+ ~ K+ (gi + g 2 )) 2 
3 " K+J+ 

( J+93 + ^+ (gi + g 2 )) 2 

U3 = 

After solving (7) with respect to / and Q as a function of 91, 92, 93, J + and x" and using the 
mass shell condition I + I~ = q\ one obtains: 



(ysx" j+ - (9-1 + 9" 2 ) 2 ) (/ + ) 2 - (9-1) 2 {J + f 

- (^7sx" J + - 2 91 • (91 + 92)) 1+ J+ = 
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which has real solutions when A^ 2 ) > 0, where 

A( 2 ) = (^x"j + ) 2 -Aq 1 -q 2 ^x"j + +4((q 1 .q 2 ) 2 -f 1 ql) 
The condition A^ 2 ) > is satisfied only if $( 2 ) > 0, where 

<j>(2) = i _ 2 (£Lgg±lgij l^ 2 D 
which gives the kinematical limits for x" and for K + 

„ _ 2{q l -q 2 + | ft \ \ q 2 |) _ J+* 



One hence obtains I and Q: 



< 1 



/+ = y/ZJ+x" -2{q l -(q 1 + q 2 )) + JW) 

2 ' ^j+ x //_(^- 1+ ^) 2 



r 
Q- 



VW) + 2 (ft ■ (ft + fc)) - J+ x" 



2 J+ 



j+ v / A( 2 T + 2 (g2-(gi + g2))-V^J + a // 
+ ^sJ+x" -2 (q 2 ■ (ft + g 2 )) 



2 J+ 

while the Mandelstam invariants s 2 = (I + + p") 2 , i 2 = (J — I) 2 , u 2 = (7 — Q) 2 are given by: 

s 2 = ^~sx" I + 



*2 



^2 



/+ J+ 
2 



(/+?2 + Q+ft 



/+ Q+ 

By solving (6) with respect to x and P as a function of x', I and ft one obtains: 



x'/+ 

X 



^~sx' -I- 



a 2 

p+ _ (/ J 



x' — I 

p- =^~ sx ' -r 
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The lower limit of x' is derived when requiring x to be smaller than 1: 



I' 

r > r ■ = 



~"min 



so we introduce 



$(D = 1 - 

x' (^-/+) 

while the expression of the Mandelstam invariants s\ = (h + p) 2 , ii = (I — h) 2 and ui = 
(P- /i) 2 is: 



fl = —yfsx'l 
U\ = —y/sxP' 
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Figure 1: Transverse momentum spectrum of partons produced in p + Pb collisions at y/s = 
5.5 , 8.8 ATe^and rj £ [2.4, 4], using po = 2GeV, and kf actor = 2. The solid lines are the result of 
implementing exact kinematics in the calculation, the dashed lines refer to the case of approximate 
kinematics. 
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Figure 2: Cronin Ratio in the case of exact kinematics (solid line) and in the case of approximate 
kinematics (dashed line) with the two different choices of energy scale. 
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Figure 3: Cronin ratio in d + Au — > ir° X collisions at s/s = 200 GeV . The solid line refers to the case of 
exact kinematics and the dashed line to the case of approximate kinematics. Left and right panels refer 
respectively to y = and r\ = 3.2. The data are form ref. [11] and [27]. 
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Figure 4: Centrality dependence of Cronin ratio with and without energy loss implementation. Compar- 
ison with experimental data presented by PHENIX. 



